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Abstract 

The contour integration technique applied to calculate the opti- 
cal conductivity tensor at finite temperatures in the case of inhomo- 
geneous surface layered systems within the framework of the spin- 
polarized relativistic screened Korringa-Kohn-Rostoker band struc- 
ture method is extended to arbitrary polarizations of the electric field. 
It is shown that besides the inter-band contribution, the contour in- 
tegration technique also accounts for the intra-band contribution too. 
Introducing a layer-resolved complex Kerr angle, the importance of 
the first, non-magnetic buffer layer below the ferromagnetic surface 
on the magneto-optical Kerr effect in the Co | Pt m multilayer system 
is shown. Increasing the thickness of the buffer Pt, the layer-resolved 
complex Kerr angles follow a linear dependence with respect to m only 
after nine Pt mono-layers. 
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1 Introduction 



The magneto-optical Kerr effect (MOKE) was discovered by Rev. J. Kerr 
in 1876. Nowadays, the MOKE occurring in multilayers is investigated be- 
cause of obvious technological implications for high-density magneto-optic 
recording media ([Bertero and Sinclair, 1994] ; |Mansuripur, 1995) ). 

The first ab-initio calculation of the absorptive part of the optical con- 
ductivity tensor based on the Kubo linear response theory (|Kubo, 1957| ), 
was carried out for Ni by Wans and Callaway (1974 ). The dispersive part of 
the optical conductivity tensor, needed to get the theoretical MOKE, is usu- 
ally computed using the Kramers-Kronig relations ( paalderop et ai, 1988|) . 
However, the first magneto-optical Kerr spectra, namely those for Fe and 
Ni, were calculated on the basis of a relativistic band-structure method by 
Pppeneer et at. (1992| ). 

Commonly, the MOKE calculations for multilayers are performed using 
conventional band-structure methods and super-cells ( Quo and Ebert, 1995 ; 
Uba et ai, 19961) . The contour integration technique developed by two 
of the authors ( [Szunyogh and Weinberger, 1999]) used in connection with 
the spin-polarized relativistic screened Korringa-Kohn-Rostoker (SKKR) 
method (jSzunyogh et a/., 1994j , |1995fc [Ujfalussy et ai, 19951 ), on the other 
hand, is a more realistic approach for the magneto-optical Kerr spectra cal- 
culations in case of surface layered systems. 

In the next section (Section the theoretical framework is revisited and 
extended to arbitrary polarizations of the electric field. In Section ^72] it is 
demonstrated that besides the inter-band contribution to the optical conduc- 
tivity tensor, the contour integration technique flSzunyogh and Weinberger 



|1999j ) includes also the intra-band contribution. In Section it is shown 
that the Luttinger formalism ([Luttinger, 1967|) - and accordingly the contour 
integration technique too - provides the complex conductivity tensor even 
in the zero frequency limit for a finite life-time broadening. The symmetric 
part of the dc conductivity tensor, however, cannot be obtained by contour 
integrating in the complex energy plane. Recently, the authors discussed in 
details schemes of how to control the accuracy of the computation ([Vernes 
et al., 20UD| ), therefore in Section [3] only the main aspects of these numerical 



methods are briefly reviewed. In Section ^ the layer-resolved complex Kerr 
angles are introduced. This concept is then used to study separately the im- 
pact of the ferromagnetic surface layer and of the non-magnetic buffer layers 
below the surface in the Co | Pt m surface layered system. In addition, the 
change of MOKE with the thickness of the buffer Pt-layer for a particular 
frequency is also analyzed. Finally, the main results of the present work are 
summarized in Section [5| 
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2 Theoretical framework 



When a time dependent external electric field is applied to a solid, currents 
are induced, which create internal electric fields. The total electric field is 
then the sum of all these fields and in the long-wavelength limit is written 



as ( Eykholt, 1986|) 



E(f, t) = E e 1 ^-^) e 5t 



(1) 
-oo 



the positive infinitesimal 5 implying the field to be turned on at t = 
( |Lax, 195M )- 

Because the solid is not a homogeneous system, in linear response the 
resulting total current density is given by QMahan, 199CQ 



j M (f,t) = J2j d3r ' J dt ' s 



/lvK r,r';t,t')E„(r',t' 



(2) 



Here £ Mi ,(r, r'; t, t') is the non-local, time dependent conductivity and ex- 
presses the linear response of the current at t in r and direction \i to the 
local electric field applied at t' in r' and direction v ( |Butler et ai., 1994| ). 
The Fourier transform of Eq. (0) 

J M (<f, u') = 6? t j8(u' ^^(q, u)E 0u (ji, v = x, y, z) , 

V 

introduces the wave-vector and frequency dependent complex conductivity 
tensor S(g, uS) (Eykholt, 1986 ), which can be evaluated using the well-known 
Kubo formula ( [Kubo, 1957] ) 



£^(g» = lim - / dte^ +i ^ 

8^0+ V Jq 



1 I dx (r,4(t + mx)) 

JO \ y / eq 



(3) 



Here V is the crystalline volume, (3 = (k^T)" 1 and (. . . ) eq the ensemble 
average in the equilibrium state at t = — oo, respectively. Notice that when 
the electric field is not transverse, i.e. qEo ^ ( Eykholt, 1986| ; [Mahan 



1990| ), the Kubo conductivity tensor oS) differs from that derived from 



the Maxwell equations ( |Kubo, 1966j ), i.e. 

S M (g,w) = E(g,a>) 
where q is the unit vector along q. 



1 q-'E(q,uj)-q 



2.1 Luttinger formalism 



Based on the contour deformation method ([Yu-Kuang Hu, 1993| ), it can be 
easily shown that 



dx (r^Kit + ihX) 

„ \ / e<i 



dt' 
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where it') is the Heisenberg operator: 



i 



For a diagonal representation of the one-electron Hamiltonian H and for the 
equilibrium density operator 

p eq = f(H) ' 



e P{JK-n) + 1 



with f(H) being the Fermi-Dirac distribution function, \x the chemical po- 
tential (in the following the temperature dependence of the chemical potential 
is neglected, i.e. \i equals the Fermi level e-p), Eq. @ can be rewritten as 



= Jim _ J £mn J -q,v,n 



dt e ^ +i5)t / dt' e ^ Em - £n)t ' , 



(4) 



where 



J ^mn = H J 3 \ n ) , J-Znm = H J -<?\ m ) 



and e m denotes the eigenvalue of H corresponding to the eigenvector \m). 

Both integrations in Eq. (|J) with respect to t' and t, respectively, can 
be performed taking advantage on the Laplace transform of the identity 
( |Abramowitz and Stegun, 1972| ). Proceeding by this manner, the factor 



e St introduced in Eq. ([!]) to describe the interaction of the solid with its 
surroundings, can be seen also as a simple convergence factor ([Lax, 1958Q . 



However, the complex conductivity tensor a for finite £ = fuv + i5, i.e. 

m,n 

extends further the meaning of S. Here, 5 can be seen as the finite life- 
time broadening, which accounts for all the scattering processes at finite 
temperature, which are not part of a standard band-structure calculation. 
Introducing the current-current correlation function as ( [Mahan, 1990|) 

UfivyLliS,) V Z-^i e — e _l_ f q,rnn J -q,nm ; l u J 

V m,n tm fc « i " t ' 



finally results the well-known Luttinger formula ( [Luttinger, 1967| ) 



E y (gO = g " (g0 ~ g " a0) ■ (7) 
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A rather technical point to be noticed is that in a previous paper ( jSzun- 



yogh and Weinberger, 1999|) , the variable ( has a slightly different definition, 
namely it was introduced as short-hand notation for u + i5/h. For this mean- 
ing of C, one has to use Eqs. (2) and (3) from ( [Szunyogh and Weinberger 



1999| ) instead of Eqs. (D and (0). 



The Luttinger formula (^) differs in form from that originally deduced by 
Kubo using the scalar potential description of the electric field QKubo, 19571) . 
In fact, Eq. (^j can be straightforwardly obtained starting from the vector 
potential description of the electric field ( |Callaway, 1974|) . Nevertheless, 



the formulae arising from the different descriptions of the electric field are 
completely equivalent (|Yu-Kuang Hu, 1993| ). 

The inter-band contribution to the complex conductivity tensor £™* er (g, Q, 
can be obtained directly from Eqs. (0) and (|?p, imposing for the former 
£ m 7^ E n . Considering the intra-band life-time broadening equal to that 
taken for the inter-band contribution, Eq. (||) provides in the e m — > e n limit 
immediately the intra-band contribution 

where 

< ra (?1 = - lim <v(£0), (9) 



because cr^ ra (q, Q = 0, if ( ^ 0. In contrast to the phenomenological Drude 
term, which is supposed to give the intra-band contribution (|Oppeneer ei 



ai, 1992"!) , S™* ra (g, () is added to both, diagonal (// = v) and off-diagonal 
(/i 7^ v) elements of the complex conductivity tensor ( Pppeneer, 1999| ). In 
conclusion, calculating a^ u (q,() i n accordance with Eq. (|), both the inter- 
and intra-band contribution is present in the Luttinger formula (|7[), i.e. 

C) = Ej£ er (?, C) + (q, C) . 



2.2 Contour integration technique 

The contour integration technique, originally developed for q = by [Szun- 
yogh and Weinberger (1999|) , is extended below to calculate the current- 
current correlation function for arbitrary polarizations of the electric field. 
This technique permits the evaluation of Eq. @ by performing a contour in- 
tegration in the complex energy plane at finite temperature. The integration 
can be performed by exploiting the properties of the Fermi-Dirac distribution 
in the selection of the contour flMahan, 1990p , 
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with z6C. Namely, that f(z) is analytical everywhere, except the so-called 
Matsubara poles ([Nicholson et al, 1994j) 



z k =e F + i{2k- 1)5t , k = 0, ±1, ±2, . . . 

(St = irksT) and for energies parallel to the real axis situated in-between 
two, successive Matsubara poles ([Wildberger et al, 19951 ), f( e ± 2iA$r) = 
f(e). 

Initially, two contours I\ and T 2 are considered around the eigenvalues 
e m and e n including a finite number of Matsubara poles, see Fig. g. The 
contour parts of I\ parallel to the real axis (z = e ± iSj, j = 1, 2) are taken 
in-between two, successive Matsubara poles 

(2Nj - 1)5 T < 8j < (2Nj + 1)5 T , (10) 

with Ni being the Matsubara poles in the upper and N 2 in the lower semi- 
plane, respectively, included in r 1; e.g. 8j = 2Nj5x for j = 1,2. T 2 is 
Ti mirrored along the real axis, hence includes N 2 Matsubara poles in the 
upper and N\ in the lower semi-plane. In order to not have e n — ( and e m + ( 
contained in the contour T 2 and Ti, respectively, the only constraint applied 
is 

5 > S 2 . (11) 



Using the residue theorem, it was shown that flSzunyogh and Weinberger. 



Sm-Sn + C, 2tt J Tl (z-e m )(z-e n + () 



5 T A 1 



+ 1 > 



and 



e m - + C 2vr (z - e n ) (z - e TO - C) 



2 

v 2 



+ z — / 



Observing that as long as ( ^ 0, the sum of these two expressions vanishes, 
when e m = e n , no restriction regarding to the eigenvalues must be imposed, 
such that by using the resolvent ( | Weinberger, 1990]) 



\n)(n\ 

J 

Z Sri. 
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immediately follows that 



<v(?> = i> dz fi z ) ^(?; 2 + C> *) - p dz f( z ) &iw(q) z,z-C) 

Jt 1 Jr 2 



-2i& 



r 



' JVi N 2 

d-fiu{q; z k + (, z k ) + 2J z k, z k - 

k=-N 2 +l fc=-iV 1 +l 



Here due to the trace, the quantity 

~ /- h 
<V(g; zi, z 2 ) = ~ ~ Tr 



2irV 



4 G{ Zl ) r, g{z 2 ) 



'121 



(13) 



contains the inter-band (e m 7^ e n ) and the vanishing intra-band (e m = e n ) 
contribution too, i.e. 



(7 



10=0 



inter / 



■ 0- 



(14) 



Thus the contour integration technique - through Eq. (O) - preserves all 
the features of the current-current correlation function introduced by the 
Luttinger formalism (Section |2.1| ). 

Originally, ^^(q; z±, z 2 ) was introduced for q = as an auxiliary quan- 
tity to be evaluated for the calculation of the residual resistivity (u, T = 0) 
of substitutionally disordered bulk systems ( [Butler, 1985| ). Since its exten- 
sion to deal with inhomogeneously disordered layered systems ([Weinberger el 
at., 1996|), this quantity nowadays is widely used to calculate the magneto- 



transport properties of multilayers ([Blaas et at., 1999|) . As it was shown in 
a previous paper ( |Szunyogh and Weinberger, 1999j ), it also plays a central 
role, when the magneto-optical properties of semi-infinite layered systems 
are calculated. In the case of electric fields with an arbitrary polarization, 
&nv{q] z%, z 2 ) obeys the following symmetry relations 

cv(<f; z i, z 2) = &up(-qi z 2, z l) 
a^(q; z* 1} z 2 ) = a* Ufl (q; z u z 2 ) = a* v (-q; z 2 , z x ) , 



dz f(z) au„(-q; z - (* , z ) 



by which Eq. (|T2| ) can be written as 



% 



dz f(z) duv{q\ z + (,z) - 



ri 



-2i5 T ^ [ a ^(q]Z k + C,z k )+al lv (-q]z k -C,z k )] . (15) 

k= -N 2 +l 

Although <7^((f, 0) is obtained directly from Eq. (^j) by taking ( = 0, 
its evaluation requires right from the beginning a single contour integration, 
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namely that along Ti, because the residue theorem provides directly 



2ir J T (z-e m )(z-e n ) 



.07 



TV £-~< 



^From this expression results immediately the inter-band contribution cr^ er (<f, 0), 
when e m ^ e n and in the e m — > e n limit this leads to 

,0/M 1 / dz M f .St A 1 
de n 2tt J Tl (z - e n ) 2 vr fe= £^ +1 (z k - e n ) 2 

Therefore, according to Eq. (p~3|) , 

, Ni 

<V(£0) = <X> dz f(z) a^ iu (q;z,z) - 2i5 T ^ *iw(<fi %k, Zk) (16) 

k=-N 2 + l 

includes the intra-band contribution too, i.e., based on Eq. (^j), 

*M o) = < er (£ o) - < ra (<?1 • ( 17 ) 

In conclusion, it can be stated that within the contour integration tech- 
nique, the optical conductivity tensor obtained by using the Luttinger for- 
mula (|7|), contains both inter- and inter-band contributions, i.e. based on 
Eqs. (HI), ADD and (§, 



inter/ - f \ _ a intev(g> q\ 



2.3 Static and sharp bands limit 

In order to derive the dc electrical conductivity, it has been shown, that first 
the limit of q — ► has to be taken ( [Luttinger, 1964| ). Eqs. ([]) and (|7|) lead to 
finite values for q —>■ in the static regime (u — ► 0), or, separately, on sharp 
bands limit (5 — > + ). It should be recalled that the contour integration 
technique can be only used for uj = calculations, when 



*V (5) 



duo 



uj=0 



because of the constraint ([TT|) applied to the contour I\, the sharp bands 
limit (5 = 0) is unreachable even when T = K. 

Alternatively, in the static limit and for finite life-time broadening, Eq. 
(P~5| ) is given by 

^ (*) = \ [ES? (5) + S« (5)] - \ [E« (6) - EH 
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where 



The sharp bands limit of its symmetric part flLuttinger, 1967|) 

S« = ~ [Eg + E^] , with E<*> = lira Eg) (6) , 

on the other hand, exists, but is restricted to the intra-band contribution 
and is given by the Landau form of the dc conductivity QGreenwood, 19581 ) 

df(e r 



V(s) = ^ V 



de r 



As a consequence of the presence of the Dirac 5-function in this expression, 
EjS can be calculated performing an energy integration along the real axis 



when T > 0, i.e., 



v(s) 



/"OO 

/ cfe 






/ —OO 



where cr^ (e) is the arithmetic mean of the four (21,22), -21,2 = e ±i0, 
given by Eq. (TL3|) in the <f — > limit (|Butler, 1985|; |Weinberger et ai., 1996|) . 



3 Computational details 

In principle, the contour 1^ along which one integrates in Eqs. (|I5|) and flip]), 
should extend from +00 to —00. In practice, however, the lower limit of 
the contour T\ is set to the bottom valence band energy e^, i.e., no core 
states are assumed to contribute. For the upper limit, on the other hand, 
s u = s~f + £&bT (£ G N) instead of 00 is taken, because the Fermi-Dirac 
distribution decays fast. (For more details, see Fig. |l[) 

In calculating the difference between Eq. (|TJ) and flTfi|), one can distin- 
guish between four different contributions, i.e. 

4 

The part of contour Ti in the upper semi-plane contributes 

i>e u +iSi 

<T$(<lX)= / dz/fc) [&p,(qiz + C,z)-ap,($z,z)] 

Je h +iO 

- / dz f(z) a^i-q; z - (*, z) , 

Je b +iO 
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whereas the contour part in the lower semi-plane contributes 



+ 



E^—iO 
e n —iS2 



dz f(z) [cy,((f; z + (, z) ~ Ziwifi z, z)] 



dz f(z) a^ u (-q;z- (*,z) 



e\>— iO 



It should be mentioned, that in a previous paper the sign of the latter was 
misprinted, see Eq. (25) from ( |Szunyogh and Weinberger, 1999| ). The Mat- 
subara poles have two contributions: one coming from the Ni — N2 poles 
situated in the upper semi-plane 



= -2iS T ^2 [ a »M\ z k + C, z k) ~ cry»/(gS z k, z k ) 

k=N 2 +l 

+^{-q; z k - c , z k ) ] 

and an other one from the 2N2 poles near and on both sides of the real axis 

C) = -2i5 T ^2 [^(9; z k + C, z k ) - a^(q; z k , z k ) + a* v (-q; z k - (*, z k ) 
fc=i 

+°nv{q; z* k + C, z* k ) - a^(q; z* k , z* k ) + a*^(-q; z* k - (*, z* k )] ■ 



In the present paper, the current-current correlation function needed to 
evaluate cr$(q,() is computed using Eq. (I3|), relativistic current operators 
(|Weinberger et at., 1996|) and the Green functions G(z) obtained by means 



of the spin-polarized relativistic screened Korringa-Kohn-Rostoker (SKKR) 
method for layered systems QSzunyogh et at., 1994| , |1995| ; [Ujfalussy et al.\ 



1995|) . Because of the finite imaginary part of the complex energy variable z, 



the calculation scheme includes the so-called irregular solutions of the Dirac 
equation too ( [Szunyogh and Weinberger, 1999j ). 

^From the computational point of view, besides the Matsubara poles, 
the optical conductivity tensor T,^ u (q,() as given by Eq. (|7|), depends also 
on the number of complex energy points n z considered for the energy inte- 
grals involved (07^ '(q, Q), on the number of fc-points used to calculate the 
scattering path operator and the a^ u (q;z ± friu + i8, z) for a given energy 
z, respectively. Recently, the authors have proposed two schemes to control 
the accuracy of the z- and /c-integrations computing £ Miy (<f, C) QVernes et at. 



2000| ). F° r that reason, only the main aspects of these numerical methods 



are given below. 

The accuracy of the integrations with respect to z is controlled comparing 
the obtained results by means of the Konrod quadrature ( Laurie, 1997 ; Ual- 
vetti et at., 200U| ), /C 2nz+ i<j^j((f, £), with those computed by using the Gauss 
integration rule ( |Press et at., 1992] ), Q ng .cr^ £), on each contour part in 
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both semi-planes (j = 1, 2) - for the notations used, see |Vernes et al. (20001) . 
Hence along a particular contour part, cxjf2(<f, () is said to be converged, if 



max 



Kan.+*T$tf, - Gn^M I < e„, (18) 



where e z is the accuracy parameter. One advantage of this scheme is that the 
integrands have to be evaluated only in 2n z + 1 points, because the 2n z + 1 
Konrod-nodes include all the n z Gauss-nodes. Another advantage is, that 
Eq. (|T8| ) is fulfilled for any, arbitrary small e z , as test calculations performed 
for q = have shown flVernes et al, 2000| ). 

In order to compute the involved two-dimensional fc-space integrals with 
arbitrary high precision, a new, cumulative special points method was devel- 
oped by the present authors ( |Vernes et al., 2000]) . This method exploits the 



arbitrariness of the mesh origin (|Hama and Watanabe, 1992| ) and requires to 
evaluate the integrands only for the fc-points newly added to the mesh. Test 
calculations for q = have shown QVernes et al, 20001) , that a requirement 



similar to Eq. flTB]), i.e. 

max \ S ni afj, u (q; z 1 , z) - Sn^a^q; z' , z) \ < eg, (19) 

for any z on the contour or Zk Matsubara pole (z' = z + (,z — (*) can be 
fulfilled with arbitrary high accuracy eg. 

Furthermore, it was shown ( |Vernes et al., 2UUU ), that if the z— and k- 



integrations are performed and controlled in the manner presented above, 
the computed optical conductivity ^^(q, () does not depend on the form of 
the contour in the upper semi-plane. Hence our computational set-up for 
S^(g, () is completely specified by e z , eg and the number of Matsubara poles 
A^ 2 near and on both sides of the real axis. (The latter is taken in accordance 
with the life-time broadening 5, i.e., to fulfill the condition: 2<5 2 = 5.) 



4 Results and discussions 

The layered system studied in the present paper consists on a mono-layer of 
Co on the top of fcc-Pt(lOO) with Pt-layers serving as "buffer" to bulk Pt 
QPustogowa et al, 1999| ): 

Co | Pt m = C0 ( i) | Pt (2 ) I ... I Pt (n _!) | Pt (n) | Pt (bulk) 

with the subscript in parenthesis being the layer index. The bottom valence 
band energy eb was taken at —1 Ryd, e u = Ep + 8/cbT and the Fermi level £f 
is that of Pt bulk, i.e. —0.039 Ryd. (Pt bulk acts as a charge reservoir for 
the layered system.) 

The optical conductivity calculations were carried out for q = 0, T = 300 
K and using a life-time broadening of 0.048 Ryd, i.e. N2 = 2 Matsubara poles 
near and on both sides of the real axis. The computation is less influenced by 
the Matsubara poles considered in the upper semi-plane, as it was already 
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mentioned in Section |||, therefore we have taken Ni — N 2 = 35 poles to 
accelerate the computation in the upper semi-plane. The convergence criteria 



10' 



a.u. 



flT|) and (JD were fulfilled for e 2 = e g 

It was shown ( |Weinberger et ai., 1996| ), that in case of layered systems 
the cr^lq; z±, z 2 ) over which one has to integrate, see Section |2?|, can be split 
into intra- (p = q) and inter-layer (p 7^ q) contributions. This decomposition 
of the current-current correlation function in Eq. ([13]) , makes the optical 
conductivity for layered systems to be of the form 



P =i 

where the layer-resolved optical conductivities are given by 

p = 1, . 



(20) 



££,(9,C)=EWC) 

9=1 



n 



(21) 



Because in the present work, q — 0, 5 and T are fixed, in the following these 
variables are omitted and the optical conductivity tensor is simply denoted 
by E Ml/ (a;). In the case of the layered system Co | Pt m , for the inter-layer 
contribution it is verified with an accuracy of 10~ 15 fs" 1 that 

= . 

It was also found that these inter-layer contributions are always smaller than 
the intra-layer contributions TP^uj). 

When linearly polarized light is reflected from a magnetic solid, the re- 
flected light becomes elliptically polarized and its polarization plane is ro- 
tated with an angle 9k with respect to the incident light. The former effect 
is characterized by the ellipticity Ek and the phenomenon is called magneto- 
optical Kerr effect (MOKE). In the polar geometry, considered in the present 
work, both the incident light and the magnetization are perpendicular to the 
surface. For a precision up to several degrees the complex Kerr angle is given 
then by ( |Reim and Schoenes, 1990 | ) 



00 



(22) 



The frequency dependence of 9k and Ek in case of Co | Pt m surface layered 
system for different number of Pt buffer layers (m = 0, . . . , 3) is shown in 
Fig. |2|. As a general trend of these curves, it can be observed, that apart from 
the m = situation at low photon energies, our calculations are in agreement 
with the experimental fact, that the MOKE decreases, if the thickness of the 
buffer Pt increases ( |Cao et at., 1998| ). 

In fact, m = means, that one considers only the surface Co-layer in the 
calculations, i.e., Eqs. (EO) and (|2"TD are taken for n = 1. Although in this 
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case, the layer-resolved optical conductivity equals S jtM/ (c<;), the correspond- 
ing Kerr spectrum is the most peculiar one in comparison with the results 
obtained for m > 1. Once the contribution of different Pt buffer layers be- 
low the surface are accounted for (m > 1), the complex Kerr angle changes 
dramatically in the whole frequency range and the calculated MOKE for 
Co | Pt m= 3, see right panels in Fig. 0, possesses already the features known 
from the experiments: #k has two typical, local minima about 2 and 4 eV 
and £k has a flat region in-between two, local extrema in the middle of the 
frequency range ( |Uba et ai, 1996Q . The frequency width, however, is smaller 



than that known from experiments. Below 1 eV and above 5 eV our theo- 
retical spectra are richer in fine details than the experimental data ( |Ebert ei 



ai, 199S| ). Finally, it should be mentioned that, both the Kerr rotation angle 
and ellipticity are approximatively three times smaller than those measured 
for thin films of Co-Pt alloy flWeller, 1996ft . 

To understand this, we introduce layer-resolved complex Kerr angles by 
using instead of £^(0;) in Eq. (|22"D , the layer-resolved TP (uj), i.e., 



fcP _ QP , j£ p _ ^y^) 



1 + — ££> 



XX \ 



(p=l,...,7i) . (23) 



For Co | Pt m=3 , is plotted in Fig. [|. In contrast to a homogeneous 
system, like bcc-Fe or fcc-Co ( [Huhne and Ebert, 1999| ), it is found that the 
surface resolved MOKE cannot predict correctly the complex Kerr effect in 
our inhomogeneous, surface layered system Co | Pt m=3 , see left panels of Fig. 
[| Inspecting the layer-resolved complex Kerr angles of the three Pt-layers 
(right panels in Fig. |^), it can be seen that both 9^ and of the first Pt- 
layer below the surface Co-layer is in the experimental range ( [Weller, 1996|) . 



This shows that the first buffer layer is as important as the magnetic surface 
in an inhomogeneous layered system. Thus in order to exploit the relatively 
big MOKE in the first buffer layer, one has to bring Pt into the surface 
Co-layer, e.g. by alloying, and has to prepare thin films. This theoretical 



finding is completely in-line with the known experimental facts QHatwar ei 
ai, 1997| ). 



in 



In addition, we have also studied the variation of MOKE in Co | Pt 
with the thickness of the Pt buffer for hu) = 0.68 eV. The change in the 
layer-resolved MOKE for m up to 15 is given in Fig. |j. The Co | Pt m 
system has only six Pt-layers below the surface included self-consistently 
and the other Pt-layers are all bulk-like. This construction is based on the 
fact, that no changes occur neither in MOKE nor in the optical conductivity, 
if the system has more than six self-consistently included Pt-layers. (Test 
calculations were carried out for Co | Pt m= g and ftw = 0.68 eV by changing 
the number of bulk-like Pt-layers in the system.) 

As can be seen from Fig. f|, the layer-resolved complex Kerr angles up 
to nine (ten) Pt mono-layers do not depend linearly on the thickness of the 
Pt buffer. This is in agreement with other, ab-initio ( piaas et ai, 2000]) 
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and model ( [Bruno et ai, 1999f ) magneto-transport calculations (u = 0) for 
different multilayer systems, but does not confirm the situation found per- 
forming super-cell calculations ( Perlov and Ebert, 200C ), where linearity of 
the surface resolved optical conductivity seems to occur after a few layers. 

We have found, that the Kerr angle of the surface Co-layer increases 
with the thickness of Pt, whereas that of the first Pt-layer below the sur- 
face, decreases and the opposite holds for the layer-resolved Kerr ellipticity. 
Although, these two layers provide the main contributions to the MOKE in 
Co | Pt m , they alone are not sufficient to describe the complex Kerr angle of 
the whole system (over nine Pt-layers), because at least twelve Pt-layers are 
needed to get #k and Ek stabilized around -0.005 and 0.009 deg, respectively. 
The change in the surface resolved MOKE with respect to m is compensated 
mainly by the change in the layer-resolved MOKE arising from the first non- 
magnetic layer below the surface, and hence #k and Ek finally are converging, 
for simplicity the layer-resolved MOKE of the Pt-layers below the third one 
are not shown in Fig. |4j. They are situated always in-between the results 
obtained for the surface Co- and the first Pt-layer below. The deeper the 
Pt-layer, the smaller its MOKE is, such that the contribution arising from a 
Pt-layer below the twelve-th layer has a really neglectable influence on the 
total MOKE. 

Finally, in Fig. |^ we show that both, the layer-resolved and the total 
optical conductivity tensor, respectively, can indeed be computed for to = 
with a finite life-time broadening by means of the contour integration 
technique, see also Section |2.3| . As can be seen from Fig. [5], the imaginary 
parts of all TP (u = 0) and T^ v {u} = 0) vanish (with an accuracy of 10 -3 fs _1 ), 
whereas the real parts of the same quantities remain finite. However, these 
values cannot be used to calculate the MOKE based on Eqs. fl22|) and ([23|), 
because these expressions are diverging in the static limit. 



5 Summary 

The contour integration technique applied to Luttinger's formalism for the 
optical conductivity tensor has been extended to electric fields with arbitrary 
polarizations and is shown to straightforwardly account for both the inter- 
and intra-band contribution. Hence within our technique there is no need 
to approximate the intra-band contribution by the so-called semi-empirical 
Drude term. The optical conductivity tensor for finite life-time broadening 
can be computed also in the zero frequency limit, however, the dc conductiv- 
ity cannot be obtained integrating in the complex energy plane. The contour 
integration technique was also completed here in deriving formulae for the 
current-current correlation function at zero frequency and vanishing life-time 
broadening (£ = 0). 

Introducing the concept of the layer-resolved complex Kerr angles and 
computing the magneto-optical Kerr effect for the Co | Pt m surface layered 
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system in polar geometry, we have shown that (1) the layer-resolved complex 
Kerr angle for the ferromagnetic surface alone cannot predict correctly the 
magneto-optical Kerr effect in inhomogeneous multilayers. In the case of 
Co | Pt m the first Pt-layer below the surface is as important for the polar 
MOKE as the ferromagnetic surface Co-layer itself. (2) The Kerr spectra of 
Co | Pt m=3 possess all the known, experimental characteristics, but in order 
to exploit the relatively big complex Kerr angle arising from the first Pt-layer 
below the surface, one has to consider a Co-Pt alloy at the surface. 

In addition, it was also shown that the layer-resolved MOKE shows a 
non-linear dependence with respect to the thickness of the buffer Pt up to 
m = 9 for Co | Pt m ; only for m > 12 the complex Kerr angle is fully 
converged. 

Finally, we have demonstrated that the contour integration technique de- 
veloped for the Luttinger formalism can be used straightforwardly with a fi- 
nite life-time broadening in the zero frequency limit. Calculations performed 
provide a purely real optical conductivity tensor at zero frequency. 
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Figure 1: Selection of contour i\ and T 2 . e m and e n are eigenvalues of the 
Hamiltonian, ( = hu + i<5, with u being the frequency and 5 the life-time 
broadening. Ep is the Fermi level, the bottom valence band energy and 
e u = Sj? + £&bT (( 6 N). Zfc = ef + % {2k — 1) 5t are Matsubara poles with 
k = 0, ±1, ±2, . . . and 5t = nksT. 5 1)2 were selected to fulfill the condition 
given by Eq. (|10|). 
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Figure 2: Complex Kerr effect for the Co | Pt m layered system with m = 
(o), 1 (•), 2 (□), 3 (■) Pt-layers below the surface Co-layer (left panel) 
and separately for m = 3 (right panel). 5 = 0.048 Ryd, T = 300 K and 
e z = er = 10~ 3 a.u., respectively. 
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Figure 3: Layer-resolved complex Kerr effect for Co | Pt m=3 . In the left 
panel the Kerr spectra resolved for the surface Co-layer (o) are compared 
with that of the layered system (■). The right panel show the complex Kerr 
effect arising from the Pt-layers: ▲ refers to the Kerr spectra resolved for 
the first Pt-layer below the surface, * for the second Pt-layer and T for the 
last Pt-layer, respectively. 
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Figure 4: Total and layer-resolved complex Kerr angle for fko = 0.68 eV 
and the Co | Pt m multilayer (m = 0, . . . , 15), containing below the surface 
Co-layer six self-consistently included Pt-layers followed by bulk-like layers. 
(The used symbols represent the same type of data as in Fig. |3|.) 
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Figure 5: Real and imaginary part of the total and layer-resolved optical 
conductivity for Co | Pt m= 3. The optical conductivity resolved for the surface 
Co-layer is given by o and that resolved for the Pt-layers by A, * and T, 
respectively. The optical conductivity of the multilayer is represented by *. 
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